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Abstract 

We study the impact of star formation and stellar feedback prescriptions on galaxy properties predicted by means 
of "stripped-down" versions of independently developed semi-analytic models (SAMs). These include cooling, star 
formation, feedback from supernovae (SNe) and simplified prescriptions for galaxy merging, but no chemical evolution, 
disc instabilities or AGN feedback. We run these versions on identical samples of dark matter (DM) haloes extracted 
from high-resolution iV-body simulations in order to perform both statistical analysis and object-by-object comparisons. 
We compare our results with previous work based on stripped-down versions of the same SAMs including only gas 
cooling, and show that all feedback models provide coherent modifications in the distribution of baryons between the 
various gas phases. In particular, we find that the predicted hot gas fractions are considerably increased by up to a factor 
of three, while the corresponding cold gas fractions are correspondingly decreased, and a significant amount of mass is 
ejected from the DM halo. Nonetheless, we also find relevant differences in the predicted properties of model galaxies 
among the three SAMs: these deviations are more relevant at mass scales comparable to that of our own Galaxy, and 
are reduced at larger masses, confirming the varying impact of stellar feedback at different mass scales. We also check 
the effect of enhanced star formation events (i.e. starbursts modes), defined in connection with galaxy mergers. We find 
that, in general, these episodes have a limited impact in the overall star formation histories of model galaxies, even in 
massive DM halos where merger-driven star formation has often been considered very important. 
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1 Introduction 

In order to understand the complex process of galaxy for- 
mation over the entire cosmic history of the Universe, 
many different physical processes need to be taken into 
account. These processes act on different scales and their 
interplay appears critical for an appropriate description of 
the chain of events leading to the build up of present-day 
galaxy population. However, our comprehension of the 
physical processes acting on the baryonic components of 
these haloes is still limited. 

A number of theoretical methods have been introduced, 
trying to get a better understanding of galaxy forma- 
tion and evolution: among these different methods semi- 
analytic mod els (hereafter SAMs - for a review, see e.g. 
Baugh 120061) of galaxy formation and evolution have be- 
come a widely used tool, thanks to their flexibility and 
(relatively) low computational costs. In these models, rel- 
evant physical processes are included by assuming empir- 
ically and/or theoretically motivated prescriptions, cou- 
pled through a set of differential equations that describe 
the mass and energy flows between the different galac- 
tic components (i.e. halo, bulge and disc) and baryonic 
phases (i.e. stars, hot and cold gas). 

A number of competing models have been proposed, 
assuming different (but equally plausible) descriptions of 
the relevant processes. Although the analysis of discrep- 
ancies and similarities between predictions from indepen- 
dently developed SAMs has been the subject of a num- 



ber of recent studies (see e.g. iFontanot et al.ll2009l 1201 1 



De Lucia et alj 1201 ll and discussions therein), the role 
played by the overall SAM "architecture" (i.e. the more 
technical details of the construction of the models, and 
the interplay of their different assumptions) has not been 
analysed in detail. 



A fi rst step in this direction was given in lDe Lucia et al 



J20T0I tereafter DL10). In this study, we used "stripped- 
down" versions of three independently developed SAMs 
including only gas cooling and galaxy mergers. The mod- 
els were run on identical sets of merger trees extracted 
from A-body simulations, so as to remove any systematic 
effect due to the adopted description of the dark matter 
evolution. We found a reasonable level of agreement be- 
tween predictions from the three models, for the physical 
processes considered. In particular, the agreement is very 
good at dark matter haloes (DMHs) scales comparable to 



those of our own Milky-Way (MW). In the other hand, 
for larger masses, corresponding to those of clusters at 
z = 0, we found significantly different results in the pre- 
dicted amount of cold gas, largely due to the different as- 
sumptions for the hot gas distribution inside DMHs, and 
to different treatments of the "rapid cooling" regime. 

In DL10, we also compared the different treatment 
for the dynamical evolution of substructures and galaxy 
mergers, based either on fitting formulae derived from 
numerical simulations or analytic models accounting for 
dynamical friction, tidal stripping and tidal shocks. We 
showed that these different assumptions result in signifi- 
cant differences in the timings of mergers, with important 
consequences for the formation and evolution of massive 
galaxies. 

In this research note, we extend the analysis of DL10 to 
investigate the influence of different prescriptions adopted 
for the physical mechanisms of star formation and SN 
feedback in both "quiescent" and "starburst" modes, 
avoiding to consider other processes like metal evolution, 
AGN feedback and disc instabilities and using simplified 
treatment of galaxy mergings. We stress that, given our 
limited understanding of the physical processes consid- 
ered, all the prescriptions we will discuss in the following 
are equally plausible, so our analysis is not aimed at iden- 
tifying the "best" model for star formation and feedback. 
Rather, the our aim is to analyse the influence of different 
model ingredients (and of the various modelling that can 
be adopted for specific processes) on the predicted proper- 
ties of galaxies and their redshift evolution. As in DL10, 
we run our models on the same sets of merger trees ex- 
tracted from numerical simulations. Results from our pre- 
vious study allow us to control residual differences due to 
a different model for gas cooling and galaxy mergers. 

In this wor k we use stripped-down vers ions of the Mu- 
nich model bv lDe Lucia & Blaizol d2007l) . the GALACTI- 
CUS model of Benson] ( 2012 ) as an extension of Durham 
m odel of iBower et all (120061) . and the MORGANA model 
of iMonaco et al.l ( 12007b . In all models we implement ra- 
diative cooling of a gas with primordial composition, star 
formation and SN feedback in "quiescent" and "starburst" 
regimes, including models of "super-wind" ejection from 
the DM halos and later re-accretions, and simplified pre- 
scriptions for galaxy mergers. We do not modify the 
choice of parameter values with respect to the original 
calibrations: we do not then expect predictions from these 
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stripped-down versions to be in any way representative of 
real galaxies, since they miss some key physical processes 
by construction, namely AGN feedback, metal evolution 
and disc instabilities. On the other hand, some of the ex- 
cluded physical mechanisms are known to have a strong 
impact on predicted galaxy properties: since each model 
fine-tuning is done on versions including these additional 
processes, they could mask or reduce any difference due 
to the processes considered here. 



2 Models 

In the following, we briefly review the merger trees used 
and describe the main ingredients of the semi-analytic 
models focusing on the physical processes considered. 
We refer the reader to the original papers for more detailed 
descriptions; the not-interested reader may skip directly to 
sec [3] 

It is worth stressing, that we make no effort to reduce 
the differences between the "cooling only" realizations 
shown in DL10. This choice allows us to keep our predic- 
tions as close as possible to the original formulation of the 
three SAMs under investigation. However, we consider 
two different sets of predictions relative to the Durham 
model, by considering both an isothermal and a cored pro- 
file for the hot halo. This choice allows us to compare 
predictions of the fiducial Durham model with those of 
a model that gives results closer to those obtained using 
the fiducial cooling models in MORGANA and the Munich 
SAM. Moreover, in order to exclude differences due to 
merger times (t mlg ) assigned to satellite galaxies, we fo- 
cus on a "no-merger" realization assuming t mlg = oo for 
all satellites. We also consider "instantaneous merger" 
realizations (< lmg = 0), and we will comment on the pre- 
dictions of these models whenever appropriate. 

We choose the same Chabrier IMF for all models and 
we switch off metal production in our model realizations 
(i.e. primordial composition is assumed during the entire 
evolution of our model galaxies). As a check, we have 
rerun our realizations allowing metal enrichment: the re- 
sults presented in Section [3] are modified in the expected 
direction (e.g., cooling rates are systematically increased), 
but our main conclusions remain valid. 



2.1 The Simulations and Merger Trees 

In this work we take advantage of merger trees ex- 
tracted from two large, high-resolution cosmological sim- 
ulatio ns, namely the Mill ennium Simulation (MS here- 
after, Springel et al. 20051) and the Millennium-II Simu- 



lation (MSII hereafter, iBovlan-Kolchin et alj|2009l) . The 
MS follows N = 2160 3 particles of mass 8.6 x lO 8 M /h 
within a comoving box of size 500Mpc/h on a side. 
The MS-II follows the evolution of the same number of 
particles in a volume that is 125 times smaller than for 
the MS (lOOMpc/h on a side), with a correspondingly 
smaller particle mass (6.9 x 1O 6 M0 /h). For both sim- 
ulations, the cosmological model adopted is a A cold 
dark matter (CDM) with Q m = 0.25, Q b = 0.045, 
h = 0.73, n A = 0.75, n = 1, and a s = 0.9. The Hub- 
ble constant is parametrised as H = lOOhkm/s/Mpc. 
Group catalogues were constructed using the standard 
friend-of-friend (FoF) algorithm, and each group was then 



proce ssed using the algorithm SUBFIND dSpringel et al 
20011) to identify self-bound substructures. 



We then consider the FoF merger trees constructed as 
detailed in DL10, and the same two sets of trees consid- 
ered there. A first sample (the "MW-like" sample) has 
been constructed by selecting from the MS-II 100 haloes 
with logio(AI-2oo /Mq) between 11.5 and 12.5 at z = 0. 
Here, M200 is defined as the mass within a sphere of den- 
sity 200 times the critical density of the Universe at the 
corresponding redshift. A second sample of 100 haloes 
was selected from the MS by taking haloes that have a 
number density of 10~ 5 h 3 Mpc -3 at z ~ 2, and that end 
up in massive groups/clusters at z = 0. The adopted num- 
ber density has been chosen to be comparable to t hat of 



submillimetre galaxies at z ~ 2 (IChapman et all 12004 , 
see e.g.,). We thus refer to this sample as the "SCUBA- 
like" sample. 

2.2 The Munich model 

The stripped-down versio n of the Munich model use d in 
this work is built upon fhe lDe Lucia & Blaizoj (12007b im- 
plementation and uses prescriptions for the star forma- 
tion and feedback that have been described in detail in 
Croton et al.l (120061) . Cold gas is associated only to the 
disc component of model galaxies, and both a "quies- 
cent" and a "starburst" mode for star formation are con- 
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sidered. Cold gas surface densities higher than a given 
threshold £ cr it are required for quiescent star formation to 
take place. This critical value may be expressed in terms 



of the galactocentric distance (R) and the virial velocity 
of the host halo (Ki JKauffmannl 19961) : 



120 x 



Kir/200 kms~ 
R/kpc 



-M &P c- 



(1) 



which can be translated into a critical mass, assuming the 
cold gas is uniformly distributed over a disc with outer 
radius rdi SC - The disc scal e length (r s ) is c omputed using 
results from the model by iMo et al.l (119981) . and the outer 
radius of the disc is assumed to be r^- lsc = 3 x r s . The 
critical gas mass for quiescent star formation to take place 
is: 



M crit = 3.8 x 10 



9 ^vir 



r<n s 



200/fcms- 1 lOKpc 



M . (2) 



The star formation rate (SFR hereafter) ip is then assumed 
to occur at the rate: 



<P 



mun — "mun 



^cold — -frfcrit 

_mun 

'dyn,D 



(3) 



where the star formation efficiency is set to a mun = 0.07, 
and Tjy™ D = rdisc/Kir is the disc dynamical time. The 
adopted modelling leads to episodic star formation self- 
regulating to maintain a level close to that corresponding 
to the critical surface density. 

In addition to the quiescent mode, the model also al- 
lows for a "collisional starburst" mode of star forma- 



tion (ISomerville et all 120011) . triggered by galaxy merg 



ers. The amount of cold gas converted into stars through 
this mode depends on the baryonic (gas + stars) mass ratio 
of the merging objects. If it is larger than 0.3, the event 
is classified as a "major" merger, and all cold gas present 
in the two merging galaxies is converted into stars. In 
the case of a minor merger, the model assumes that only 
a fraction /b IS of all available cold gas is converted into 
stars: 



brs 



ebr 



Mi 
Ah 



(4) 



where M1/M2 represents the baryonic mass ratio (M2 > 
3Mi), and ab rs = 0.7 and &b rs = 0.56 have been cho- 
sen to provide a good fit to the numerical simulations 



of ICox et al.1 (120041) . Mergers also involve mass transfer 
between the disc and bulge components of the remnant 
galaxy. A detailed description on how galaxy mergers 



affect galax y morphology i s prese nted in iDe Lucia et al. 



( 120 1 lb and Fo ntanot et alj d201 II) . For the purposes of 
this work we do not distinguish between bulge and disc 
components of model galaxies, and focus on their global 
properties. 

As for stellar feedback, the Munich model links the 
amount of cold gas reheated by SNe (AM™°) in a given 
time interval to the mass of stars formed in the same time- 
step (AM s f): 



AMZ n 



e rht AM sf , 



(5) 



where e r ht = 3.5. The energy released in the same time 
interval (A_Z?sn) can be written as: 



A£ SN = 0.51/ s V/s N AM sf , 



(6) 



where 0.5V^ N represents the mean energy in SNe ejecta 
per unit mass of stars formed (Vsn = 630fcms _1 based 
on standard SNe theory and a Chabrier IMF), and 7?sn = 
0.35 parametrises its efficiency in reheating the disc cold 
gas. 

Adding the reheated gas to the hot halo without chang- 
ing its specific energy leads to the total thermal energy 
change (AE hot ): 



AE] 



hot 



0.5K? N xAMri7\ 



(7) 



It is then possible to define an excess energy E CKC = 
A^sn — A^hot- If ^cxc < 0, all reheated gas is con- 
fined within the halo, otherwise a fraction AM"?™ of hot 
gas mass (A/hot) is ejected from the parent halo through a 
"super-wind": 
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V SN 
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vir 



e rht AM s f • 



(8) 

where E\ wt = 0.51/^.Mhot represents the total ther- 
mal energy of the hot gas. The ejected material may be 
reincorporated at later times as the p arent DMH grows 
dPe Lucia. Kauffmann & Whitell2004l) : 



M mun = 



n M e j e 



Tdyn,H 



(9) 
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where ri™f n = 0.5 is a free parameter which controls 
the amount of reincorporation per halo dynamical time, 
7"dyn,H = -Rvir/Kir (-Rvir being the virial radius of the 
parent halo). 

2.3 GALACTICUS 

In the SAM comparison we presente d in DL10 w e made 
use of a stripped-down version of the iBower et al.l (120061) 
implementation of the Durham model. There we tested 
two versions of the Durham model, the difference lying 
in the assumptions made for the profile of the hot gas dis- 
tribution: w e defined a "standa rd" model with a /3-profile 
(as used by Bower et al.ll2006l) and a "modified" version 
using an isothermal profile. 

In this work, we take ad vantage of the new GALACTI- 
CUS code ( BensonL 2012 ) to recreate the stripped-down 
versions of the Durham models of DL10. GALACTICUS 
is designed to be highly modular to facilitate the explo- 
ration of different descriptions of key physical ingredi- 
ents. We set-up GALACTICUS to run with the same as- 
sumption^ regarding gas cooling and gas profiles as in 
the standard and modified versions of the Durham code 
used in DL10. In the following, we will refer to these 
realizations as GALACTICUS-CP (for the f3 or cored pro- 
file) and GALACTICUS-IP (for the isothermal profile). We 
explicitly check that the "cooling only" stripped down 
versions of GALACTICUS reproduce with good approx- 
imation the predictions of the "cooling only" Durham 
model presented in DL10. We then include the same 
treatment for star formatio n an d stellar feedback as de- 
pict ed in ICole et all d2000h and Benson et alJ (120031 see 
also Bower et al ,l2006hT in older to create an equivalent of 
the Durham model needed for our analysis. Note that both 
the GALACTICUS-CP and GALACTICUS-IP models differ 
significantly from the default model in the GALACTICUS 
toolkit. 

Galaxy sizes play a key role in this model since they 
determine dynamical times and rotation speeds in discs 
and spheroids. Given the angular momentum of cooling 
gas, the radii of galactic discs are computed by solving 



In the GALACTICUS realization used in this work, all other rele- 
vant assumptions, including the definition of formation times, DM pro- 
files and concentrations, merger time calculation, galaxy size calcula- 
tions, and major merger definitions are treated as in the original Durham 
model. 



for the equilibrium radius at which rotation supports them 
against gravity in the combined potential of disc, spheroid 
and NFW dark matter halo (including the effects of adi- 
abatic contraction). When spheroids are formed through 
major mergers (see below) their radii are computed by as- 
suming conservati on of (internal plus orbital) energy. Full 
details are given in lCole et al. I fcOOOl) . 

Star formation is assumed to occur in galactic discs, 
their cold gas reservoirs being depleted at a rate deter- 
mined by the star formation timescale T s d f ur : 



¥>dur = A , i C oid/T" s d f Ur . 



(10) 



The star formation timescale is a function of the disc 
circular velocity U d ™ taken at the half-mass radius : 
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dur dyn,D 



/ disc 



V 200fcms- 



/3dur 



(ID 



whereT dyu,D 



r HR 

disc 



/^disc re P resen ts the dynamical time 
of the disc, while a^ur = 0.0029 and ^ur = —1-5 are 
free parameters. 

When two galaxies merge, the merger is deemed to 
be "major" if the less massive of the galaxies is at least 
30% of the mass of the more massive galaxy. In such 
cases, the stars of both galaxies are redistributed into a 
spheroid. Major mergers always trigger a starburst. Mi- 
nor mergers may also trigger a starburst if the gas fraction 
in the more massive galaxy exceeds 10%. In starbursts, 
the gas content of the merging galaxies is placed into the 
spheroid component of the merger remnant, where it pro- 
ceeds to form stars on a timescale given by eq. [TTlbut 
with adur = 0.5, /?dur = 0, and with a dynamical time 
and velocity computed for the spheroid, using the method 
described bv lCole et al. 

Stellar feedback is modelled by assuming that the rate 
at which cold gas is reheated (M^ t r ) is directly related to 
the SFR: 



Vdur 



lit 



V. 



HR 



Vdi: 



(12) 



where Vhot = 485fcms~ 1 and 7 = 3.2 are free parame- 
ters. The reheated gas is not instantaneously returned to 
the hot phase, but it is stored into a separated reservoir, 
M TSV . This material is then added back to the hot phase at 
a rate equal to: 
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dur 
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dur Jw rsv 



(13) 



r dyn,H 

where f?™}" 1 = 1.26 and Td yni H is the dynamical time of 
the dark matter halo. It is worth noting that the reheated 
material in GALACTICUS behaves as the ejected material 
in the Munich model, i.e. as a separated gas phase, which 
does not take part in the exchange of energy and mass in- 
side the parent DMH, and which is slowly reincorporated 
into the hot phase. For a sake of simplicity, in the fol- 
lowing, we will then refer to this component as ejected 
material: 



M* 



M, 



rht 



(14) 



2.4 MORGANA 



The MOdel for the Rise of GAlaxies aNd Agns ( MOR- 
GANA) was first presented in iMonaco et al.l (120071) . The 
treatment of star formation and stellar feedback in this 
code fol lows the results of the multiphase model for the 
ISM by Monaco! ( 2004 ). For the purposes of this work 
we consider the combination of parameter values adopted 



when using a Chabrier IMF (ILo Faro et al.L 120091) . 



Discs are treated as "thin" systems: SN remnants blow 
out of the disc soon after they form and most of the SNe 
energy is injected into the halo hot gas, and only a few 
percent of SNe energy injected into the ISM. The star for- 
mation timescale is of the form: 



'sf.D 



9.1 



/ Scold, D 

V im qP c- 2 



-0.73 



fcold.D 
0.1 



0.45 



Gyr, 



(15) 

where / co id,D represents the cold gas fraction (in the disc 
component). Gas re-heated by stellar feedback is ejected 
from the disc to the an external reservoir of baryons (i.e. 
the "halo" component) still bound to the host DM halo, 
at a "hot wind" rate Mhw assumed to be equal to the star 
formation rate: 



Mi 



hw,D 



(16) 



A fraction /th.D = 0.32 of SN energy E$n is carried away 
with this ejected material, so the contribution Eh w of the 
hot wind to the thermal energy of the hot halo gas is: 



Ehw,D = /th,D-E: 



SN" 



1mor,D 

msN 



(17) 



where msN represents the mass of newly formed stars per 
SN. A further energy contribution is added by assuming 
that one type la SN per year explodes each 10 12 Mq of 
stellar mass. This rather crude implementation of energy 
from SNe la does not influence much model results. 

Cold gas flows into the bulge component either by 
mergers or by disc instabilities. Moreover, a fraction of 
the cooling flow is allowed to fall directly into the bulge. 
Star formation in bulges is assumed to take place in a 
"thick" regime, where SN energy is effectively trapped 
within the ISM. In this case, SFR is assumed to follow a 
Schmidt-Kennicutt relation, with a timescale: 



( s c 



sf,B 



-0.4 



Gyr . 



(18) 



\ MqPC 2 / 

The size of the starburst, necessary to compute the gas 
surface density, is estimated as follows. Gas is assumed 
to have no angular momentum, and thus to be supported 
by turbulen ce generated by SNe . Under very simple as- 



sumptions dLo Faro et aU 12009b . the velocity dispersion 



in this case can be written as: 



Scold = CT 



Gyr 



-1/3 



(19) 



where oo = 60 km s _1 is treated as a free parameter. It is 
then assumed that the size of the starburst region is such 
that cr cc ,id equ ates the rotation curv e of the bulge, assumed 
to be flat (see lLo Faro et al ]|2009j for details). 

The rate at which hot gas is ejected from a bulge to the 
host halo is assumed to be: 



Mi 



hw.B — 



V hot 





Vmor.B if Vb < Vhot (20) 

ifVs > Vhot 



This is done to take into account the ability of the poten- 
tial well of a massive bulge to keep hot gas confined; the 
parameter Vhot is chosen to be 300 km/s, corresponding 
to the typical thermal velocity of a ~ 10 7 K hot phase. 
The corresponding energy carried by the hot wind is: 



E, 
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\Ahot 
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(21) 
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where we use /th,B = 0.1 for the fraction of SN energy 
carried away. 

A bulge also ejects cold gas into the halo, due to the 
same SN-driven turbulence that sets the starburst size 
(Eq. [19}. This "cold wind" from the bulge to the halo 
component is assumed to occur at a rate: 



M, 



J d.B-Pub^ub 

Rb 



(22) 



where Rb is the half-mass radius of the bulge, while P u b 
and v u b represent the probability that a cold cloud is un- 
bound (i.e. it has a velocity larger than the escape ve- 
locity of the bulge Vb) and the average velocity of un- 
bound clouds (both probabilities are computed assuming 
a Maxwellian distribution of velocities with rms a C o\d)- 

The hot and cold halo gas components keep track of 
both the mass and the energy received respectively from 
hot winds from discs and bulges and from cold winds 
from bulges. The halo receives winds both from the cen- 
tral galaxy and from satellites. Whenever the gas phases 
of the halo component are too energetic to be bound to 
the DM halo, they are allowed to escape to the IGM as a 
galaxy "superwind". In particular, if the energy £hot,H °f 
the halo hot gas mass Mhot.H overtakes the virial energy 
E v i T by more than a factor /wind = 2, a superwind occurs 
at a rate: 



M e 7 e OT H = I 1 - 



/wind-^vir 
-Ehot.H 



C,M\ 



hot.H 



R\ 



(23) 



where c s represents the sound velocity in the halo. A sim- 
ilar formula determines the ejection of cold superwinds: 



/winded; 



disp 



<t h M, 



cold,H 



(24) 



A fraction (/back = 0.5) of the mass ejected by the 
DM halo is later re-incorporated (i.e. added to infalling 
IGM), when the parent DM halo reaches an escape veloc- 
ity larger than the (thermal or kinetic) velocity the gas had 
when it was ejected. 

2.5 Comparison between the models 

All three models considered relate the SFR in the disc 
component to the amount of cold gas there available, but 



make different assumptions for the timescale of conver- 
sion. In the Munich model, the star formation proceeds 
on a timescale directly proportional to the dynamical time 
of the disc, while in Durham-like models the timescale 
of star formation is rescaled with some power of the cir- 
cular velocity of the disc. Finally, MORGANA assumes 
a star formation timescale consistent with the Schmidt 
law. The Munich model explicitly accounts for a critical 
mass threshold for star formation, while MORGANA and 
Durham-like models do not. 

Besides a quiescent mode of star formation, both the 
Munich and GALACTICUS-CP models assume an en- 
hanced star formation regime occurring during galaxy 
mergers, with a fraction of the total cold gas available be- 
ing turned into stars in a very short timescale (that of the 
model integration). In MORGANA, cold gas is associated 
also to the bulge component. Because of the complete 
loss of angular momentum, gas in bulges is concentrated 
to very small sizes. Then the higher gas surface densi- 
ties force the Schmidt-Kennicutt relation to give much 
higher SFRs and shorter star formation timescales than 
those corresponding to star formation in discs. In our ref- 
erence runs, we explicitly exclude merger events and the 
enhanced SFR by setting t mTg = oo. We will discuss 
the impact of the "starburst" mode on the predicted star 
formation histories by considering the predictions of the 
instantaneous-merger runs in section |3~2| 

In MORGANA the amount of disc gas reheated via stel- 
lar feedback (M ln t) equals the SFR (ip), while in the Mu- 
nich model these two quantities are proportional, via the 
parameter e rn t > 1. This implies that the amount of re- 
heated gas per unit star formation is always larger in the 
Munich model than in MORGANA . A slightly different 
choice has been made in the GALACTICUS models, where 
the amount of reheated gas remains proportional to ip, but 
is additionally assumed to scale as a power of the disc 
velocity. 

An important quantity in the balance of the baryonic 
content of each DMH is given by the amount of baryons 
ejected. In both Durham-like and Munich models, the 
fraction of ejected mass is directly proportional to the 
mass of stars formed in the same time interval, while in 
MORGANA the dependency of the ejection rate on stellar 
feedback is mediated via the estimate of the thermal (ki- 
netic) energy of the hot (cold) gas halo phase. All mod- 
els assume that this ejected material is reincorporated into 
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the halo at later times, with different assumptions for the 
reincorporation rates. In Durham-like and Munich mod- 
els, this is modelled as a continuous process following the 
DMH growth, while in MORGANA only half of the mate- 
rial connected to each ejection event is re-acquired instan- 
taneously when the DMH has grown enough to overcome 
its escape velocity. 

Durham-like and Munich models employ a very simi- 
lar scheme for modeling the ejected component. In both 
cases some of the feedback reheated material is excluded 
from the mass/energy flows between the cold and hot 
gas phases: but while this is only a fraction in the Mu- 
nich model (the remaining being instantaneously added to 
the hot phase), all reheated material is considered in the 
Durham-like models. The same reincorporation scheme 
is assumed, but with a faster rate in the Durham-like mod- 
els than in the Munich model. 

In the following, we choose to avoid the complications 
arising form disc instabilities^ by switching them off. 

3 Results 

3.1 Test cases 

We start by considering the evolution of the different bary- 
onic components in a few test cases. For consistency with 
DL10, we focus on the same 4 representative DMHs (their 
fig. 1 and 4). These have been chosen among the 100 
MW-like and the 100 scuba-like (two for each sample) as 
representative of a quiet mass accretion histories and of a 
large number of merger events. In the following, when- 
ever we refer to "central galaxy" for our DMHs, we will 
refer to the central galaxy of the main progenitor at the 
corresponding redshift. All other galaxies will be defined 
as "satellites". For the sake of simplicity, in this section 
we will show results only for the 2 haloes with quiet mass 
accretion histories. The other two haloes give consistent 
results, and are shown for completeness in appendix lAl 

2 In particular, disc instabilities have no direct effect on the SFR pre- 
dicted by the Munich model, since only enough stellar mass is removed 
from the disc to the bulge to restore stability. In the Durham model, at 
each instability event, the whole disc is destroyed and all its baryons are 
given to the spheroidal component, with any cold gas present assumed 
to undergo a starburst. In MORGANA a different choice has been made, 
by moving a well defined fraction of disc stars and cold gas to the bulge, 
where it forms stars on a T s ™g' timescale. 



Figures Q] and [2] show the redshift evolution of the 
different baryonic components in the haloes considered: 
from top to bottom, we show the cold gas fraction asso- 
ciated with the central galaxy, the hot gas fraction in the 
halo, the net cooling rate (ttincr; see below), the stellar 
mass of the central galaxy and the mass ejected from the 
DMH. In the three upper panels, we compare the predic- 
tions of the models including star formation and feedback 
(darker colours - solid lines) with previous results for the 
"only cooling" realizations considered in DL10 (lighter 
colours - long dashed lines). 

We define mNCR as the net rate at which cold gas and 
stars are deposited into the galaxy: 

• _ Af C oid(*2) + AW(t 2 ) M cold (ti)+A/ Btar (ti) 

TOn cR - t^T, ti=tl ' 

where t\ and are the cosmic epochs corresponding to 
two consecutive snapshots. The net cooling rate differs 
from the intrinsic cooling rate, since it takes into account 
both the effect of feedback in removing part of the cold 
gas from the system and the effect of star formation in 
locking some material in long lived stars. For the stripped- 
down SAM versions considered in DL10 these two quan- 
tities coincide by construction. 

The inclusion of feedback from star formation has the 
net effect of reducing the amount of cold gas available 
in MW-haloes (figure Q] this holds also for instantaneous 
merger realizations) and increasing the hot gas fraction 
with respect to the "cooling only" runs. In the MW haloes, 
this increase is particularly relevant in the Durham-like 
models, that predicts higher hot gas fraction with respect 
to the other two models. This shows that the feedback 
scheme implemented in the Durham-like models is the 
most efficient (among those considered here) in reheating 
the cold gas in the haloes at these mass scales. The cold 
gas fraction associated with the central galaxy is charac- 
terized by a marked decline at lower redshift, and predic- 
tions from the three models considered are more differ- 
ent than in the "cooling only" realizations. In our models 
there are two competing effects able to deplete the cold 
gas reservoir: the formation of long lived stars and cold 
gas removal by stellar feedback. In order to disentangle 
these two effects we consider the evolution of the overall 
cold gas plus stars and still find a decrease with respect 
to the predictions of DL10 for the cold gas component. 
These findings show that the inclusion of a strong stel- 
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MORGANA Munich GALACTICUS-CP GALACTICUS-IP 




Log(l+z) 



Figure 1 : Redshift evolutio n of the baryonic con tent for a MW-like representative DMH (with quiet mass accretion 
history; fig.l right panels in De Lucia et al. 2010b . Upper row. cold gas fractions associated with the central galaxy; 
second row. hot gas fractions; third row. net cooling rates. Blue, red, yellow and green lines refer to MORGANA , 
GALACTICUS-CP , GALACTICUS-IP and the Munich model. Dark solid line s refe r to the models considered in this 

201(lh 



work, while dashed lines refer to the "cooling only" models (IDe Lucia et al 



dashed lines) and ejected masses (red solid lines). In all models we assume t E 



Lower row. stellar masses (blue 
oo (see text for more details). 
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Figure 2: Redshift evolutio n of the baryonic cont ent for a SCUBA-like representative DMH (with quiet mass accretion 
history; fig. 4 right panels in lDe Lucia et al J2010b . Symbols, line styles, and colours have the same meaning as in fig. 03 
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lar feedback is the main driver of the different evolution 
of the baryonic components for the models considered in 
this study. 

The combined effect of the decreased cold gas and in- 
creased hot component is shown in the third row. In gen- 
eral, the evolution of rfiNCR is n °t smooth and character- 
ized by epochs of negative contribution, i.e. time intervals 
dominated by outflows able to reduce the cold gas content 
of the central object. For the MW-like haloes, the evolu- 
tion of 7TIncr show large deviations from the "cooling- 
only" configuration, with extended redshift ranges char- 
acterized by low or negative rates. This is particularly 
evident for the Durham-like realizations, while the Mu- 
nich model gives predictions that are closest to the cooling 
rates obtained in the "cooling only" runs. 

In the lower panels of fig. Q] we show the redshift evo- 
lution of the baryonic mass ejected from the DMH (red 
solid line) and the stellar mass associated with the cen- 
tral galaxy (blue dashed line). The first quantity is cru- 
cial to keep track of the total amount of baryons in each 
main progenitor, and gives important hints on the impact 
of stellar feedback. At this mass scale, the ejected compo- 
nent dominates over the stellar mass depo sited in the cen 



tral g alaxy for MW-like haloes, (see also iDe Lucia et al 



20041) . We stress that the ejected mass in the Durham-like 
realizations is comparable to the results for the other two 
models for the MW-like haloes, confirming that the differ- 
ences in the cold and hot gas fractions are mainly due to 
the different feedback efficiency and not to an enhanced 
fraction of material excluded from the mass/energy flows. 
Indeed, the similar amounts of ejected mass in Durham- 
like and Munich models in the whole redshift range, de- 
spite rjf^ 1 3> ?7™" n , confirms that the stellar feedback 
scheme adopted in Durham-like models removes a larger 
fraction of gas from the baryonic budget of these haloes 
with respect to the Munich scheme. 

In the same panel we also show the evolution of the 
central galaxy stellar mass: at this mass scale and with 
the inclusion of only star formation and stellar feedback, 
model results are very different. The Munich model 
predicts the largest stellar masses at all redshifts, MOR- 
GANA gives slightly lower values for the stellar mass, 
and both Durham-like runs predict stellar masses almost 
one order of magnitude smaller than the other two mod- 
els. The GALACTICUS-IP model predicts slightly larger 
stellar masses with respect to the Durham-like models. 



These results for the Durham-like models are consistent 
with our analysis of the net cooling rate, and are due to a 
feedback-driven starvation of the cold gas reservoir of the 
central galaxy. We checked that these conclusions also 
hold for the redshift evolution of the mean stellar mass in 
the central galaxy (averaged over the 100 realizations in 
each DMH sample). 

Similar conclusions are reached when considering the 
representative SCUBA halo (figure |2), but with some sig- 
nificant differences. For these haloes, both the increase 
of hot gas fraction and the decrease of the cold gas frac- 
tion with respect to the "cooling only" versions are less 
marked, and the prediction of different models are some- 
what closer. This is due to the fact that stellar feedback is 
more efficient in affecting the thermal state of the gas and 
galaxy evolution on a MW-like scale than in more mas- 
sive haloes. Consistently, in SCUBA-like haloes toncr 
follows much more closely the evolution of cooling rates 
in DL10, and the outflow dominated periods are less fre- 
quent or completely absent. As already noticed in DL10, 
net cooling rates can takes up values of several hundreds 
of M Q yr _1 with the exception of the GALACTICUS-CP 
model, where spikes are much less pronounced. 

At these mass scales, the ejected component does not 
dom inate over the stellar m ass for the SCUBA-like sam- 
ple dDe Lucia et all |2004 at z < 2, and the differ- 



ent feedback models predict a rather different evolution: 
in the Durham-like runs the amount of material in the 
ejected component stays constant below z ~ 3, a more 
noisy evolution is seen for MORGANA , while in the Mu- 
nich model the ejected mass decreases rapidly and is neg- 
ligible for the baryonic budget of the haloes at present. 
The differences in the predictions of the three models are 
due to the different interplay between the ejected and rein- 
corporated fraction, and again shows that in Durham-like 
models the fraction of ejected material is larger than in 
the Munich model, leading to an almost constant ejected 
mass, despite the faster reincorporation rate. 

Finally, the stellar masses predicted by the three SAMs 
for the SCUBA-like sample are much closer than in the 
MW-like sample, with only GALACTICUS-IP predicting 
slightly higher stellar masses. We interpret also this result 
as an effect of a less efficient stellar feedback in SCUBA- 
like haloes for Durham-like models: in these models gas 
cooling is always dominating over cold gas removal (third 
row) and more material for star formation is available. 
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Therefore, in this case the marked decrease of the cold 
gas fraction associated with the central galaxy is the ef- 
fect of a more efficient star formation. 

3.2 Star Formation Histories 

We now turn our attention to the stellar mass assembly of 
the central galaxy as predicted by our SAMs. In figure [3] 
we consider the evolution, averaged over the whole sam- 
ple of 100 halos, of the stellar mass in the central galax- 
ies between two contiguous snapshots: in our reference 
models, with no galaxy mergers, this corresponds to the 
mean SFR in the object between the two snapshots (dark 
lines). In order to quantify the relative contribution of 
the "starburst" mode to central galaxy mass assembly, we 
also consider the predictions of the instantaneous merger 
run (t mrg = 0) and proceed as follows. We first com- 
pute the mean stellar mass variation for central galaxies 
in the runs with instantaneous mergers: these quantities 
include both the SFRs in central galaxies and the contri- 
bution from satellite mergers. In the no-merger runs, we 
then consider the mean stellar mass content in satellites 
already accreted by the main progenitor at each redshift. 
We then subtract this from the mean stellar mass varia- 
tions for central galaxies in instantaneous merger runs. 
This quantity (shown as light lines in the figure) does not 
strictly correspond to the SFR in the central galaxy in the 
instantaneous merger runs, since enhanced SFR episodes 
in satellites are not properly subtracted. Clearly, the in- 
stantaneous merger runs correspond to an extreme case, 
since in a realistic SAM run only a fraction of satellites is 
allowed to merge onto their central galaxy. Nonetheless, 
the comparison of dark and light line provides a conserva- 
tive upper limit to the overall starburst mode contribution 
to the assembly of the central object. 

In each panel of figure [3] we consider the mean varia- 
tion of stellar mass in the two different sets of runs, aver- 
aged over the 100 MW-like (upper panels) and SCUBA- 
like haloes (lower panels). In the left panels, we di- 
rectly compare the no merger runs with the instantaneous 
merger runs modified as described above for each SAM, 
to show the contribution of the starburst mode in each con- 
figuration. In the right panels, we compare the different 
predictions of our SAMs in the same runs. 

Let us focus first on the MW-like sample. As far as 
the no merger runs are considered, MORGANA and the 



Munich model predictions are very close to each other at 
early times, and start diverging at lower redshifts, with 
MORGANA predicting less star formation than the Munich 
model. At the same mass scales, the GALACTICUS real- 
izations show lower SFRs at all redshifts. With respect 
to the Munich model, the difference is about one order 
of magnitude, while at lower redshifts the difference with 
MORGANA predictions is reduced. The GALACTICUS-IP 
predicts systematically higher SFRs at all redshift than the 
GALACTICUS-CP run. These results are consistent with 
our findings in sec. 13.11 and are directly related to the dif- 
ferences in net cooling rates among the SAMs. We com- 
pare these results with the analogous prediction for the 
instantaneous merger runs, thus showing the maximum 
expected contribution from the starburst mode of star for- 
mation. In general the overall growth rate of the central 
galaxy is enhanced in all SAMs. The largest contributions 
are seen at low redshift for MORGANA (at z < 3, reaching 
a factor of 8 at z ~ 0) and for GALACTICUS-IP (a factor 
of 4 at z < 1), while the Munich model shows a smaller 
(at most a factor of 2) increase at all redshift. The smallest 
modifications are seen for GALACTICUS-CP (only a few 
percent). 



We then focus the SCUBA-like sample. For the no 
merger runs, the predictions of the different models are 
generally closer than for the MW-like sample. At these 
mass scales, the predictions of the different GALACTI- 
CUS realizations are more similar to those from the Mu- 
nich model and MORGANA , and they show higher SFRs 
at earlier times (z > 2). Again, GALACTICUS-IP realiza- 
tions predicts systematically higher SFRs with respect to 
GALACTICUS-CP realizations. Also for this sample, the 
instantaneous merger runs predict an enhancement of the 
stellar mass formed, but in this case the variations are well 
below a factor of two in all cases, with the maximum devi- 
ation seen for MORGANA at intermediate (z ~ 3) redshifts 
(it is worth recalling that the two examples of figures|2]and 
|7]show spiky deposition rates for three out of four models; 
these differences are not visible in the SFR when it is av- 
eraged over 100 halos). These results shows that, at these 
mass scales, the impact of merger-driven starburst is lim- 
it ed. This is in l ine w ith the phenomenological estimate 
of lSargent et al.l (12012b of an 8 - 14 per cent contribution of 
mergers to the total SFR. 
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Figure 3: Average evolution of central galaxy stellar mass as a function of redshift (over 100 DMHs - see text for 
more details). Blue, red, yellow and green lines refer to MORGANA , GALACTICUS-CP , GALACTICUS-IP and the 
Munich model. Dark dot-dashed lines refer to the i mrg = oo realizations, while light solid lines correspond to the 
instantaneous merger runs. Upper panels refer to the MW-like sample, while lower panels to the SCUBA-like sample. 
Left panels compare the no merger and the instantaneous merger runs for each SAM separately, while right panels 
compare the average prediction of different SAMs in the same configuration. 
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Figure 4: Comparison of central galaxy stellar masses for the 100 MW-like (left panel) and SCUBA-like(n'g/zf panel) 
haloes. In each panel, blue circles, yellow asterisks, green crosses, purple squares and red triangles refer to model 
predictions at z = [0; 0.5; 1; 2; 3] respectively. 



MORGANA Munich GALACTICUS-IP MORGANA Munich GALACTICUS-IP 




Figure 5: Same as fig.|4]for cold gas masses associated with the central galaxy in the 100 MW-like (left panel) and 
SCUBA-like(n'g/!f panel) haloes. Symbols and colours are the same as in fig. S] 
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3.3 Stellar and gas content at different red- 
shifts 

In order to get more insight into the mass assembly pro- 
cess predicted by our models, we compare in fig. |4]the 
model stellar masses for the central galaxy at different 
redshift on an object-by-object basis. Left and right pan- 
els refer to the MW-like and SCUBA-like haloes, respec- 
tively. Each pair of models show some degree of corre- 
lation in the predicted stellar masses, and this confirms 
the overall consistency of the models. For the MW-like 
haloes, the stellar content predicted by the Durham-like 
realizations is offset low with respect to predictions from 
the other models. Moreover, the overall slope of the re- 
lation looks steeper than the one-to-one correlation, with 
smaller galaxies deviating more from the one-to-one rela- 
tion. This shows that the difference is due to the mass 
depencence of the efficiency of feedback in regulating 
star formation. For the SCUBA-like haloes, the stellar 
masses predicted from the three models considered are 
very close. GALACTICUS-CP is the only model deviating 
considerably from the one-to-one relation, particularly at 
low redshift, reflecting the fact that these galaxies grow 
more slowly in the GALACTICUS-CP realizations than in 
MORGANA or the Munich model. We also consider the 
corresponding predictions for the "instantaneous merger" 
run (not shown in the figure): we find that all central 
galaxies in SCUBA-like haloes at all redshift lie along the 
one-to-one line, while our conclusions are qualitatively 
unchanged in MW-like haloes. For these runs, the scatter 
in the correlations is reduced in both samples at all red- 
shift: we thus conclude that the presence of mergers helps 
the predictions of the models to converge to a common 
value, and this effect gets stronger with increasing num- 
ber of mergers. 

We also consider the corresponding cold gas content 
associated with the central galaxy on an object-by-object 
basis (fig. 01. The differences between the SAMs are even 
more evident for this physical quantity in both samples. 
MORGANA and the Munich model show the strongest cor- 
relation both for the MW-like and SCUBA-like haloes, 
while the Durham-like realizations show no correlation 
at all with the predictions of the other two models. It 
is worth stressing that in SCUBA-like haloes Durham- 
like realizations predict cold gas amounts 1 — 2 orders 
of magnitude lower than the corresponding MORGANA 



and Munich predictions. The results are similar between 
GALACTICUS-CP and GALACTICUS-IP , implying that 
this depletion is not only related to the smaller cooling 
flows associated with the cored-profile, but that stellar 
feedback is also playing an important role in removing 
cold gas content from the haloes. However, it is not possi- 
ble to indicate stellar feedback as the only responsible for 
the smaller stellar masses obtained in MW-like haloes. In- 
deed, the cold gas content of these haloes in the Durham- 
like runs may be even larger than the corresponding MOR- 
GANA and Munich predictions. Therefore, it is the inter- 
play between star formation and feedback that is respon- 
sible for the different predictions. 



4 Discussion & Conclusions 

We compared predictions of three independently devel- 
oped semi-analytic models of galaxy formation, focusing 
our analysis on their assumed modelling for the physical 
processes involving star formation and stellar feedback. 
Following the same approach as in lDe Lucia et al.l (120 KX 
DL10) we define "stripped down" SAM versions includ- 
ing cooling, star formation, feedback from supernovae 
(SNe) and simplified prescriptions for galaxy merging, 
and we run them on the same samples of DMH merger 
trees, extracted from the Millennium and Millennium-II 
Simulations. 

The choice of "stripped-down" versions of the SAMs 
has the advantage of avoiding complications due to other 
processes like disc instabilities, metal enrichment and 
AGN feedback. In our "stripped-down" versions we ei- 
ther assume i mrg = oo or i mrg = 0, in order to remove 
the additional degeneracies due to different definitions of 
merging times (see DL10). Our aim is to discuss the influ- 
ence of specific model ingredients on the physical prop- 
erties of model galaxies. In the following, we summarise 
and discuss our main findings. 



• Supernovae feedback: As expected, switching on 
star formation and stellar feedback has important 
consequences on the different gas phases in DMHs, 
with respect to the "cooling only" SAM versions: 
the amount of cold gas available is reduced, while 
the hot gas fraction is increased. While we find an 
encouraging level of consistency between the mod- 
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els, the specific star formation and SN feedback pre- 
scriptions in each model induce differences in the 
stellar and gaseous content of galaxies. This is in 
line with what is found in numerical simulations 
of MW-like galaxies perf ormed with different codes 
(IScannapieco et al.l 120121) . In particular, we find that 
the scheme adopted by the Durham-like models pro- 
vides larger hot gas fractions with respect to the other 
two models and a rapid depletion of the cold gas 
fraction in galaxies. The amount of gas that is in an 
'ejected' component (i.e. temporarily not associated 
with the halo/galaxy) is similar in the three models 
considered, though its redshift evolution can differ 
significantly. 

• Stellar content: If we consider the average SFRs 
in the MW-like sample, MORGANA and the Munich 
model predictions show a good level of agreement 
at z >, 2: predictions from the two models deviate 
at later times, which has only a limited effect on the 
predicted z — stellar masses. For this sample, the 
most striking differences are between predictions of 
MORGANA (or the Munich model) and those from 
the Durham-like realizations. In fact, the star for- 
mat ion and feedback sc hemes implemented follow- 
ing Bower et al. d2006l) predict significantly lower 
SFR levels at all redshifts. As a consequence, the 
Durham-like models predicts systematically lower 
stellar masses for the corresponding central objects 
with respect to both MORGANA and the Munich 
model. The discrepancy is a weak function of stellar 
mass itself, being smaller for more massive galax- 
ies (within a given DMH sample). For SCUBA-like 
haloes, the differences between the predicted SFRs 
are smaller than for MW-like haloes, and the pre- 
dicted stellar masses are much closer at all redshifts, 
the biggest difference being a more rapid assembly 
of stellar mass in the central object at earlier times in 
the Durham-like runs. 

• Quiescent and starburst modes of star formation: 
We analysed the impact of the "starburst" mode of 
star formation (usually associated with galaxy merg- 
ers) by comparing runs obtained with infinite and 
vanishing galaxy merging times; this approach al- 
lows us to give an upper limit on the contribution 
of the starburst mode. In most cases we found an 



increase of the average mass assembly rate of the 
central galaxy not larger than a factor of 2, with 
the largest contributions found in the MORGANA 
model. This shows that the contribution of the star- 
burst mode of star formation in the overall SFR bud- 



get is limited, in line with, e.g. JSargent et al.l (120121) . 
We also find that the inclusion of merger driven star- 
bursts decrease the scatter in the predicted stellar 
masses at given redshift between haloes in the same 
DMH sample. 



These results extend and deepen our conclusion pre- 
sented in DL10 to include the comparison of different 
approaches to the modeling of star formation and feed- 
back in SAMs. Despite the general coherent picture for 
the impact of the energy injection connected to stellar 
feedback on the distribution of baryons into the differ- 
ent gas phases, each feedback model is characterized by 
its unique pattern in tracing the redshift evolution of the 
different baryonic components (cold gas, hot gas, ejected 
gas) and this information is fundamental to understand the 
building up of galaxy properties in a cosmological con- 
text. 
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For consistency with DL10, we show in this appendix 
the two DMHs chosen among the 100 scuba-like and the 
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cant number of merger events along their assembly his- 
tory (compare with results in sec. I3.lt . 
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Figure 6: Redshift evolution of the b aryonic content of a M W-like representative DMH (with merger-dominated mass 
accretion history; fig. 1 left panels in iDe Lucia et al Symbols, line styles, and colours are as in fig.Q] 
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Figure 7: Redshift evolution of the baryo nic content of a SCU BA-like representative DMH (with merger-dominated 
mass accretion history; fig.4 left panels in lDe Lucia et al 1 l2010l) . Symbols, line styles, and colours are as in fig.Q] 
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